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The three-dimensional real scalar model, in which the Z2 symmetry spon- 
taneously breaks, is renormalized in a nonperturbative manner based on the 
Tamm-Dancoff truncation of the Fock space. A critical line is calculated by 
diagonalizing the Hamiltonian regularized with basis functions. The marginal 



^ ■ {(j)^) coupling dependence of the critical line is weak. In the broken phase the 



canonical Hamiltonian is tachyonic, so the field is shifted as (pix) <p{x) -\-v. 

The shifted value v is determined as a function of running mass and coupling 

so that the mass of the ground state vanishes. 
PACS numger(s):11.10.G, ll.lO.H, ll.lO.S, 11.30.Q 



Typeset using REVT[TiX 



*e-mail : taka@higgs.phys.kyushu-u.ac.jp 
''^e-mail : yahiro@fish-u.ac.jp 



I. INTRODUCTION 



Although the field theory has been studied for a long time, we have no complete un- 
derstanding of relativistic bound state problems. We have to solve the quantum chromo- 
dynamics (QCD) nonperturbatively to investigate low-energy hadronic physics. The most 
reasonable nonperturbative method for such a problem is the lattice gauge theory. In the 
method, it is easy to calculate the lightest particle state, but not to evaluate the excited and 
scattering states. Although there exist many nonperturbative prescriptions in the Hamilto- 
nian formalism, the formalism has been abandoned so far since Lorentz invariance and the 
renormalizability are not obvious. 

The light-front (LF) Tamm-Dancoff field theory [|l| is known as a method based on 
the Hamiltonian formalism. The Hamiltonian is constructed by quantizing fields on the 
light-cone and truncating the Fock space. The truncation is called the Tamm-Dancoff ap- 
proximation 0. The completeness of the Fock space is approximated as J2n=i ~ 1) 
where A^td ~^ oo corresponds to the full theory. It is possible to calculate mass spectra 
nonperturbatively by diagonalizing the Hamiltonian in the space. The resulting mass spec- 
trum seems to be accurate for low-lying states, since the pair creation and annihilation are 
suppressed in the LF field theory 0, especially the pair creation of particles from the LF 
vacuum is kinematically prohibited. The prohibition warrants the truncation to be reliable. 
This is an advantage of this method, but, at the same time, it causes a problem. No pair 
creation means that the vacuum is trivial H\0) = 0, i.e., the vacuum is always symmetric in 
the LF field theory. This causes serious problems, because various phenomena are explained 
as results of spontaneous symmetry breaking (SSB) in the relativistic field theory. How 
can the LF field theory describe SSB? It has been said that the zero-mode operator leads 
the LF field theory to SSB 0,0. We have a constraint equation for the zero-mode, if the 
system is quantized in a box with a periodic boundary condition. New interactions induced 
by the zero- mode affect the renormalization P|, independently of whether the symmetry 
breaks spontaneously or not. A relation between the zero- mode scenario and the conven- 
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tional understanding based on the effective potential in the equal-time field theory is not 
clarified yet. There are several works |^ for the two-dimensional real scalar model, in which 
the zero-mode constraint equation is solved under reasonable approximations to describe 
SSB of the Z2 symmetry. At least for the two-dimensional case, the light-front formulation 
seems to be consistent with the equal-time one, if the zero-mode contributions are treated 
properly. This approach to SSB has never been applied to higher dimensional models having 
ultraviolet divergences. 

According to our previous work based on perturbation, the canonical Hamiltonian 
yields a tachyonic spectrum for the broken phase owing to renormalization effects. We have 
to expand the field around the correct vacuum expectation value (VEV) to restore positivity 
of the spectrum. There exists an RG invariant relation among VEV and running parameters, 
and the VEV is determined from the relation as a function of the original parameters. The 
broken phase is an RG invariant surface containing the critical line in the parameter space. 
The symmetry breaking interactions were originally introduced by Perry and Wilson in the 
coupling coherence study PJTI1|], but the instability of the canonical Hamiltonian was not 



mentioned. 

In this paper, we treat the real scalar model in three dimensions and propose a nonpertur- 
bative RG prescription based on the Tamm-Dancoff truncation. The theory is regularized by 
expanding wavefunctions in terms of basis functions in a truncated momentum space. The 
mass spectrum is then obtained by diagonalizing the Hamiltonian within a space spanned by 
the basis functions. In actual calculations we take the Fock space up to three-body states. 
It is discussed through wave functions whether the truncated space is large enough. The 
marginal coupling dependence of the critical line is also analyzed. A proper Hamiltonian for 
the broken phase is constructed by shifting the field as (^[x) —>■ ip{x) -|- f . A relation of v to 
the original parameters are nonperturbatively determined by searching the critical surface 
which gives a vanishing mass for the ground state. 

The paper is organized as follows. In section ||, a qualitative aspect of mass spectra 
in the three-dimensional real scalar model is shown with perturbative Bloch-Horowitz RG 



analysis. In section |T|, a regularization based on basis functions is introduced. In section pV| , 
numerical RG is carried out by using the nonperturbative regularization mentioned above. 
The section |V| is devoted to discussions. 

Notational conventions are summarized as = {x^ ±x^)/y/2, time x~^, space 

{x~,x±), momentum k = {k~^,k±), metric = = —g^^ = —gii = 1 and others 

= 0, and LF energy = e(k) = {kj_ + r)/2k'^, where r is a mass parameter. 

II. SHIFTED FIELD AND POSITIVE DEFINITE HAMILTONIAN 

Let us consider a Hamiltonian regularized with some cut-off A„ in order for a divergent 
theory to be well-defined. The Einstein-Schrodinger (ES) equation for the Hamiltonian is 

i/(r„,A„,...;A„)|^) = E|^), (2.1) 

where r„, A„, . . . are running parameters belonging to the model. We need a certain trans- 
formed Hamiltonian which has a lowered cutoff A„4.i(A„+i < A„) and gives the same eigen- 
value E, 

i/'(r„+i,A„+i,...;A„+i)l^') =^l*')- (2-2) 

We are now looking for a prescription to get transformed parameters r„+i, A„+i which include 
radiative corrections. To do that, we use the following projection operators O], 



V = eiAl^,-Ml,), Q = e{Ml,-Al^,), (2.3) 

where Mint is a total invariant mass (Mj^^. = P^) of an intermediate Fock state. The 
transformed Hamiltonian, which gives a correct spectrum in the V space, can be written 
formally as, 

H' = VHV + VVQ—^—-QVV, =V\^). (2.4) 

E — H 

This Hamiltonian H' is called the Bloch-Horowitz effective Hamiltonian and used in many- 
body problem for the purpose of reduction of the degrees of freedom which are too much 



to treat it directly |T2[. We expect that this Hamiltonian effectively describes physics near 



the ground state. Expanding the term 1/{E — H) with respect to the coupling constant, 
integrating over the higher invariant mass space Q and rescaling the lowered cutoff A„+i to 



the original one A„, we can write down RG equations for all parameters |T3[. The effective 
Hamiltonian, however, contains a problem. The resultant RG equations seems to depend 
on the energy E at a sight. In our previous work we have shown perturbatively for 
the four-dimensional 0"^ model that the energy dependence of RG equation is negligible at 
one- loop order if the cutoff is sufficiently large. It's not clear whether it is true or not beyond 
perturbation. In this paper, we use the effective Hamiltonian ( p.4|) only in order to explain 
qualitative effects of RG on mass spectra perturbatively. 

Next, we derive a perturbative RG equation for mass parameter in the three-dimensional 
real scalar model. We assume that the Lagrangian density is 

^ = Id.^d'^^ - '-4>' - _ ^^6_ (2.5) 



The field is quantized as 



{x), (l){y)]x+=y+ = -^^(x" - y')6{xi_ - y±), (2.6) 



if the zero-mode is neglected. The factor 1/4 on the right-hand side of Eq. (|2.6|) is from the 
Poisson bracket 

{$(x), $(y) = i-dl + dl)5ix- - y-)5{x^ - yx), (2.7) 

where $(x) = 7r(x) — 9_0(x) is a primary constraint and tt is conjugate to 0. The canonical 
Hamiltonian is given by. 

The normal ordering is taken here. This is equivalent to renormalization of tad-pole diagrams 
by redefining parameters. We have to include various kinds of relevant and marginal interac- 
tions to renormalize all new interactions produced by RG transformations. Our Hamiltonian 



is assumed to have three interactions included in Eq. ( |2.8|) . The mass renormahzation is 
most important, because r is most relevant to mass spectra. It is the kinetic term which 
mainly contributes to mass spectra especially in the perturbative region. The RG equation 
for mass r is given with perturbation as 

B 



Tn+l = ( r„ + 



1927r' 



(2.9) 



where 



B = 



3 W2Zp. 
li=l '^i 



52(p-f:k,) fi?-E6(k, 

■t=l \ 1=1 



L = 



Ar 



A 



(2.10) 



n+l 



It is understood implicitly that the loop integral is made in the Q space. The factor L is 
produced by rescaling the new cutoff An+i to the original cutoff A„. By changing variables 
to Jacobi coordinates, 

ki = {xiK+, XiK±_ + Si), Yh=i Xi = 1, ELi ■Si = 0, 

< Xi < 1, —OO < Sj < oo, 

where all intermediate states have a common total momentum K, the integral B becomes 

dxidx2dsids2 1 



B 



(2.11) 



XiX2{l -Xi- X2) M2 - Ml^{xi, X2, Si, S2) ' 

where M is an external mass (an eigenvalue of Hamiltonian) and Mint a mass in the inter- 
mediate Fock state, 



Ml^^{xi, X2, Si, S2) 



gj +^« ^ sl + rn ^ (si + 52)^ + r„ 



(2.12) 



Xi X2 1 — Xi — X2 

It is assumed that the external mass M is small compared to the cutoff scale. The intermedi- 
ate mass Mint is large because it sits in the higher mass space Q. Then, M^ is much smaller 
than Mj^^.. This assumption is reasonable for the purpose to calculate a critical line and 



draw a phase diagram. We can observe from B being negative that lim„^oo = cxd if the 
initial value tq is sufficiently large and lim„^oo = ~oo if the initial value tq is small. That 
is, there exist a critical line r = rc(A) in the first quadrant of the two dimensional parameter 



space (r, A), and r > rc(A) and r < rc(A) correspond to symmetric and asymmetric phases 
respectively. In the asymmetric phase, the canonical Hamiltonian ( |2.8| ) is tachyonic, i.e., 
square of eigenmass is not bounded from below although Hamiltonian should be positive 
definite. Positivity of the spectrum is restored by expanding the field (j) around the correct 
VEV V. Substituting (j){x) = (p{x) + v into Eq. (|2.8|), we have 

Ha = Ho + Hu (2.13) 
Ho = J d'x^({d^^y + g2v'), (2.14) 

H,= /rf^xE V' (2.15) 
J ~n n\ 



i=3 



where the new mass parameter is 



g2^r + — + ^, (2.16) 



and the new coupling constants are 



fi-s = (A + -— )t;, g4 = X + —-, g5 = wv, ge = w. (2.17) 

D 2 

The V is just a free parameter at this stage. We have three parameters in the symmetric 
Hamiltonian Hs and four in the asymmetric one Ha- The RG parameter space for Ha is 
larger than that for Hs- The v should be a function of the original three parameters belonging 
to Hs, that is, there should exist an RG invariant relation among the four parameters. 
Actually, following our previous work for the real scalar model in four dimensions, we 
can find out the same RG invariant relation 6r + Af ^ = also for the three-dimensional 
case, with perturbative renormalization at one-loop order, if the 0^ interaction is neglected. 
The neglect of the marginal operator will be justified by numerical analyses shown later. 
The RG invariant relation forms a surface in the parameter space of if a- All asymmetric 
Hamiltonians on the surface describe physics for the broken phase. All running parameters 
on the surface go to zero in the limit A ^ which corresponds to infinite iterations of RG 
transformation, satisfying the RG invariant relation 6r + Af ^ = 0. The Hamiltonian Ha 
thus gives massless spectrum in the limit. All Hamiltonians on the surface converge to the 



same massless Hamiltonian in the limit, indicating that Ha gives massless spectrum at any 
point on the surface. It is known for the real scalar models in three and four dimensions 



that the effective potential is convex and has flat bottom in the broken phase |T^. We may 
choose any point in the flat region as a stable vacuum. It is likely that this fact supports the 
existence of a massless particle in the broken phase. We will continue our analysis with the 
fact in mind, that is, we flnd out the broken phase nonperturbatively by searching massless 
eigenvalues. 



III. BASIS FUNCTION REGULARIZATION 

The flrst principle of RG is to flnd a flow which gives the same physics. It is not 
easy to calculate the effective interaction part of the Bloch-Horowitz Hamiltonian directly 
for our practical purpose of nonperturbative renormalization and we have to avoid energy 
dependences of renormalization, then we consider another approach. In our framework, we 
can draw RG flows by calculating spectra. We regularize the LF Hamiltonian with the basis 
function regularization and try to calculate the critical line and the critical surface. 

In general, we can express an arbitrary state in the Fock space such as 

n n 

6\P - J2 k.)^n(ki, k2, . . . , k„) n a^km, (3-1) 

i=l i=0 

where P is the total momentum of the state and the wavefunction tpn is symmetric under 
exchanges of arbitrary two momenta. The limit A^td oo corresponds to the full theory 
and the state is normalized as 

(vl/(P)|v[/(Q)) = <52(p_Q). (3.2) 

We will set a certain small number to A^td, which is called the Tamm-Dancoff approximation. 
In this paper, the Fock space is truncated up to three body states (A^td = 3). According to 
the variational principle, the mass spectrum in the truncated Fock space is given by solving 
(diagonalizing) the ES equation 



"TP 1 p 

l^(P)) = E-^/ 

n=i V n! J 



i=0 
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p2|^(P)) = M2|^(P)). (3.3) 

We introduce a cutoff A and show the present regularization scheme later. All quantities 
belonging to the model are measured in units of the A and attached tilde. In units of the 
cutoff, eigenmass and running parameters are 

M = M/A, d = ^;/A^/^ f = r/A^ A = A/A, w = w, (3.4) 

then 

~gn = ^7„/A(6-")/2_ (3 5) 

The 0^ interaction is marginal as shown in Eq. ( |3.5| ). We will check the marginal coupling 
dependence of the critical line in the next section. The momentum is also redefined as 

X, = {xi, Xi) = {kt/P+, ki±/A), (3.6) 

where < Xj < 1 since kf' > 0. In this rescaled notation, the normalization condition Eq. 
( p.2|) of the wavefunction becomes 

n 

d^P - ^x,)|^„(xi,X2, . . . ,xO|' = 1, (3.7) 

i=l 

where 

P = (l,Px) = (l,P^/A), (3.8) 
and the rescaled wavefunctions are 

n~l 

?/^„(xi,X2, . . . ,x„) = (^P+A) ' V^„(ki,k2,...,k„). (3.9) 

The ES equation in the rescaled notation is given in appendix 0. It can be seen that the 
ES equation does not depend on the longitudinal momentum and the cutoff A. 

Regularization of transverse momenta and discretization of Fock space are closely related 
to each other. The regularization of the Hamiltonian is realized as a boundary condition of 
the wavefunctions. That is, we solve this eigenvalue problem, keeping the constraint that 

8 



7Vtd 

E 

n=l 



.1=1 



the transverse component of ip3 is zero at edges [Xi = —1, 1). The ES equation is regularized 
with the naive transverse cutoff, 

- A < ki± < A ^ -1 < Xi < 1, (3.10) 

since the uhraviolet divergence is caused by the transverse loop integral. The Fock space 
is discretized by expanding the wavefunction in terms of basis functions in the momentum 
space (—1 < Xj < 1, < a;i < 1), 

7VL(n)7VT(n) 



V'„(Xi,X2, . . . ,X„) = J2 «kl '5l,2,...,n 



An) 
''kl 

k=0 1=0 U=l 

where 



(3.11) 



^x, = 1, J2Xi = P±, (3.12) 

i=l 1=1 

and 5i,2,...,n is a symmetrizer. The Ni^{n) and N'Y[n) are taken sufficiently large so that the 
eigenvalue of the ground state can converge. Longitudinal basis functions are 

/fc(a;) = < x < 1, (3.13) 

for < f3{n) < 1 and transverse basis functions are 

Fi{X) = {l-X^)X^, -1<X<1. (3.14) 

The variational parameter l3{n) is tuned so that the ground state takes minimum eigenvalue. 
The behavior of the three-body wavefunction near the edges {x = 0, 1) is important because 
of the kinetic term proper to LF [0,^ . The coefficients determined by diagonalizing 

the Hamiltonian in the momentum space spanned by the basis functions. Mathematically, 
Eq. ( p. 3D has exact eigenvalues only in the case that the functional space is expanded in 
terms of the complete set of basis functions. We can get only upper bounds of eigenvalues, 
since the wavefunction is expanded in terms of incomplete basis functions in this calculation. 
However, we expect that the spectrum is described accurately for low-lying states, since 
shapes of the wavefunctions are simple; for example, the calculated wavefunction of the 
ground state has no node. 
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In this paper, a phase diagram which includes rescahng effects is evaluated, in accor- 



dance with Wilson's RG prescription [|T3|. It is not needed to vary the cutoff because the 
transformed cutoff is rescaled to the original one in this RG. Parameter sets which give the 
same eigenvalue are calculated for a fixed cutoff. 

The wavefunction renormalization is neglected in this calculation, i.e., Z^p = \. The phase 
diagram is dominated by the Gaussian fixed point, since the field is rescaled according to 
the canonical dimension. In order to find a non-trivial fixed point, we have to consider the 
wavefunction renormalization. These two phase diagrams, which have trivial and non-trivial 
fixed points, describe different theories. Our result makes sense as the theory dominated by 
the Gaussian fixed point. 

In this calculation, the total transverse momentum of the eigenstate is set to zero, = 0, 
and iVL,T(^) = 3 in Eq. ( p. Ill) , which is sufficiently large so as to give a convergent spectrum. 



IV. NUMERICAL RESULTS 

Mass spectrum is plotted in Fig. |1] for the case A = 10 and w = without introducing 
V. No multi-bosonic bound state is found in this calculation. We can observe a massless 
point fc(A = 10). The point is called the critical point. We can draw the critical line by 
connecting the critical points which are calculated for various A. The critical line is plotted 
in Fig. 1^ for five values of marginal coupling constants, w = 0,1, 10, 10^, 10"^. Values of 
two parameters f and A which form the critical lines are tabulated in table H for each w. 
The dependence of the line on w is weak. This is consistent with our intuition based on 
perturbative analysis, i.e., the dependence of the line on higher power operators in ip may 
be week. After this, we will switch off the marginal coupling {w = 0), and draw the phase 
diagram of three running parameters, v, r and A near the v = plane. 

The lowest eigenmass is zero on the line, corresponding to an infinite correlation length 
in the statistical theoretical language. The spectrum of the canonical Hamiltonian is not 
bounded from below in the left region of the critical line. We then conclude that excited 



10 



states are lighter than the ground state in the region if we solve the canonical Hamiltonian 
as it is. This is a result of the mass renormalization effect explained in the section The 
true Hamiltonian should be positive definite even in the left region of the critical line. It is 
natural to introduce the VEV in the broken phase in order to construct a positive definite 
Hamiltonian. The symmetry has to be spontaneously broken for the spectrum to keep 
positive. Figure §(a) shows a relation Xv'^ + 6r = among three parameters v, r and A, which 
is given by perturbative RG calculations at one-loop order; the detail of the calculations is 
shown in |^ for the case of the four- dimensional scalar model. We can see a critical line 
on the positive A axis, which is trivial because no renormalization effect is included. In 
Fig. 0(b), a critical surface is drawn. The surface is calculated by searching points which 
give massless spectrum in the three-dimensional parameter space {v,r, A). Positivity of the 
Hamiltonian is restored by introducing v. Compared to the surface (a) in Fig. ^ the surface 
(b) slants to the positive f direction if v is small and to the negative direction if v is large, 
as a result of renormalization effects. 

Components of the ground state wavefunction on the critical line shown in Fig. |^ are 
plotted in Fig. ^in order to confirm whether the TD approximation works well or not. One- 
and three-body wavefunction probabilities are plotted as a function of A on the critical line. 
The three-body component is very small even for comparatively strong coupling A. This 
means that the spectrum does not change near the critical line even if higher Fock space 
contributions are included in the calculation. This is true also in region far from the critical 
line, when the VEV is zero. The three-body component of the ground state wavefunction 
tends to decrease, as mass parameter increases with coupling constant fixed, because taking 
large mass parameter is effectively the same as taking weak coupling constant. 

Next, we investigate TD dependence of the critical surface. Figure]^ shows wavefunction 
components of the ground state on the intersection between the critical surface and the 
A = 50 plane. The coupling constant is set to the largest value in this calculation. One- and 
two-body components are dominant and three-body component is small also in this case. 
The two-body component increases rapidly near -y = 0.1, whereas the three-body component 
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increases slowly and keeps values around a few percent. It is expected near the v — plane 
that the Tamm-Dancoff approximation is good. 

V. SUMMARY AND DISCUSSIONS 

We have shown how to renormalize the Hamiltonian nonperturbatively for the three- 
dimensional real scalar model and calculated the critical line and the critical surface. There 
exist a region where the canonical Hamiltonian is unstable independently of perturbation. 
The mass spectrum is tachyonic in the left region of the critical hne. We have then introduced 
VEV to restore instability of the Hamiltonian and calculated the critical surface. The 
asymmetric Hamiltonian Ha. on the surface describes physics for the broken phase, if the 
phase has a massless mode. This is precisely true for SSB of any continuous symmetry, and 
the present approach for finding the broken phase is applicable for the case. Similarly to 
the present bosonic case, there seems to exist a region where Hamiltonian is unstable if we 
renormalize a mass term m.'ipip in fermionic theories. It is reasonable to understand that the 
instability is restored by introducing VEV of a bosonic field. 

The marginal coupling dependence of the critical line is weak. This result links to the 
fact that the Tamm-Dancoff approximation up to three-body state works well. The ES 
equation for ipi has no marginal interaction, but the ES equation for ip^ has that. The 
ground state has marginal coupling dependence only through the three-body wavefunction 
■03 which is extremely small as a result of nonperturbative calculation. This is the reason 
why the marginal coupling dependence of the ground state is weak. 

This RG program seems to be independent of energy eigenvalues. In order to check this, 
we need plural number of bound states. For example, we can use spin singlet and degenerate 
triplet states in a fermionic theory such as Yukawa model. It is possible to draw RG flows 
by using a spectrum of one-body fermionic state in a similar way which has been done in 
this paper. We can say that the RG is energy independent if the spectra of the spin singlet 
and triplet states are renormalized with the same RG flow. 
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APPENDIX A: OSCILLATOR EXPANSION OF THE FIELD 

The field (y? is expanded at a;"*" = with oscillators as 

dk 



roo roc 
/O J-oo J2k+{27T}^ 



where 



a(k),a^(k') 



APPENDIX B: EINSTEIN SCHRODINGER EQUATION 



The ES equation ( |3.3| ) becomes 



M2?/'„(Xi,X2, . . . ,x„ 



+ 



+ 



+ 



+ 
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n x2 

4! (2^)^^y§ 



?/;„(xi,X2, . . . ,x„) 

52(x,, + x, -E-=i 



xV'n(yi,y2,Xi, . . . ,Xj, . . . ,X 
^^(X, + Xj + Xfc-Ei=iyO 

6! (27r)4yy ^XiXjXfc 



.J , . . . , X^ J 



96 30 



X V'nlyi, y2, ya, xi, . . . , Xj, . . . , Xj, . . . , Xfc, • • • , x,,) 



in + l)\ /^5^.-E-=iy^ 



5'3 3 

3! V2(27r) V ri! ^ 

xVi„+i(yi,y2,Xi, . . . ,Xj, . . . ,x„) 

^5 10 ' 



5! V2(27r) 



(n+1)! /■^f(_x. + x,-ELiy. 



n! 



x^/'„+i(yi,y2,y3,xi, . . . ,Xi, . . . ,Xj, . . . ,x„) 
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+ ^3 6 J {n-1)\ I 5\^, + X, - y) 



3!V2(27r)V n\ Jyf^. ^ 



i<j V •"'"'J 

xV'n-l(y) Xi, . . . , Xj, . . . , Xj, . . . , X„) 



^ ^5 30 J (n-l)! ^ ^ (52(xi + X,- + x^ - y) 



5! V2(27r)3 V n! 



i<j<k 



1 "-^ J J 



x?/'„_i(yi,y2,xi, . . . ,Xi, . . . ,Xj-, . . . ,Xfc, . . . ,x„) 



4!(27r)2V n\ Jyfr{ 

xV^n+2(yi, y2, ys, Xi, . . . , X^, . . . , X„) 
+ ^6 15 ^^/ (^^T2)! fj^5^, + X,- - Eti yO 



x^n+2(yi, y2, ya, y4, xi, . . . , x^, . . . , x^-, . . . , x„) 



^ J2_ / (n-2)! r " (5^(x, + x,- + x^ - y) 
4!(27r)2V n! iy.!^, v'x^ 

X'0n-2(y, Xi, . . . , Xj, . . . , Xj, . . . , Xfe, . . . , X„) 



+ 



g6 90 / (n-2)! 1 " (5^(x, + x,- + Xfc + x;-EiiyO 
6!(27r)4V n! h iJ^^Ki ^Jx^jx^i 



xV'n-2(yi, y2, Xi, . . . , Xj, . . . , Xj, . . . , Xfe, . . . , X;, . . . , X„) 



^ ^5 ^5 ^/ (n + 3)! I ^ 52(x, - EtiyO 



+ 



5!2v^(27r)3V iy^ 

xV'n+3(yi, y2, ys, y4, xi, . . . , x^, . . . , x„) 



^5 60 /(n + 3)! 



^ 52(xi + Xj+Xfe + Xi-y) 



5! x/2(27r)3 V 

X V'„_3(y, xi, . . . , Xi, . . . , Xj, . . . , Xfc, • • • , X;, . . . , x„) 

+ ^6 3 ^^/ (^r+4)! ^ ^ 52(x, - Eti Yi) 



6!2(27r)4V n! iy^ 

xV5n+4(yi, y2, ys, y4, ys, xi, . . . , x^, . . . , x„) 



+ 



^6 180 /(n-4)! 



^ ^^(Xj + Xj + Xfc + X; + x^ - y) 



6! {2nY V n' •'y i<j<k<l<m ^JXiXjXkXlXm 

'^'4^n—4:iy J Xi, . . . , Xj, . . . , Xj, . . . , Xfe, . . . , X;, . . . , X^yi) ■ ■ ■ ) X^j), 



where e(x) = {X^ + g2)/x and 



dyj_ 
\.f=i VVi. 



n 



(B2) 



14 



this shows integrations over n-pieces of y, and the x's with the check symbol are removed 
from arguments of the wavef unctions. On the right hand side of Eq. (PT|), interaction terms 
in which the number of checked x's is larger than the number of particles n are switched off. 
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FIGURES 
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FIG. 1. Mass spectrum for the lowest eigenstate with A = 10. The spectrum is plotted as 
a function of the input mass parameter r. is a critical point, which gives the massless ground 
state. 
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respectively. The critical line little depends on the marginal coupling constant. 
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FIG. 3. (a)A surface for the relation XiP' + 6f = is figured. (b)A critical surface, which gives 
the massless ground state, is calculated by diagonalizing the light-front Hamiltonian in the Fock 
space truncated up to three-body states. 
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FIG. 4. Wavefunction components of the ground state on the critical hne; one- and three-body 
components are shown by diamonds and crosses as a function of A, respectively. Lines are intended 
to guide the eyes. 
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FIG. 5. Wavefunction components of the ground state on the critical hne; one-, two- and 
three-body components on an intersection between the critical surface and the A = 50 plane are 
shown by diamonds, crosses and squares as a function of respectively. Lines are intended to 
guide the eyes. 
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TABLES 

TABLE I. Values of mass parameter f which give massless modes on the v = plane are 

tabulated as a function of A in the cases w = 0, w = 1, w = 10, w = 10^ and w = 10^. The 
marginal coupling dependence of the critical line is weak. 



w\\ 


0.0 


10.0 


20.0 


30.0 


40.0 


50.0 


0.0 


0.000 


0.044 


0.105 


0.170 


0.236 


0.303 


10° 


0.000 


0.044 


0.105 


0.170 


0.236 


0.303 


IQi 


0.000 


0.044 


0.105 


0.170 


0.236 


0.303 


10^ 


0.000 


0.044 


0.105 


0.170 


0.235 


0.302 


10^ 


0.000 


0.041 


0.100 


0.164 


0.230 


0.296 
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